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ABSTRACT 

In an earlier study of the spiral galaxy M33, we photometrically identified arcs or outer spiral arms of intermediate 
age (0.6 Gyr — 2 Gyr) carbon stars precisely at the commencement of the Hi-warp. Stars in the arcs were unresolved, 
but were likely thermally-pulsing asymptotic giant branch carbon stars. Here we present Keck I spectroscopy of seven 
intrinsically bright and red target stars in the outer, northern arc in M33. The target stars have estimated visual 
magnitudes as faint as V~ 25 m . Absorption bands of CN are seen in all seven spectra reported here, confirming their 
carbon star status. In addition, we present Keck II spectra of a small area 0.5 degree away from the centre of M33; 
the target stars there are also identified as carbon stars. We also study the non-stellar PAH dust morphology of M33 
secured using IRAC on board the Spitzer Space Telescope. The Spitzer 8/j.rn image attests to a change of spiral phase 
at the start of the HI warp. The Keck spectra confirm that carbon stars may safely be identified on the basis of their 
red J— K s colours in the outer, low metallicity disk of M33. We propose that the enhanced number of carbon stars in 
the outer arms are an indicator of recent star formation, fueled by gas accretion from the Hi-warp reservoir. 

Key words, galaxies: evolution - galaxies: spiral - galaxies: individual (M33 = NGC598) 



1. Introduction 

Like the two dominant spiral galaxies in the Local Group 
(MW and M31), M33 is known to have a prominent warp. 
This warp is spectacular in the HI-21cm line: at many 
points the line of sight intersects the disk twice. In M33 
the warping commences at a radius of ~ 5 kpc and at a ra- 
dius of ~ 10 kpc the rotation axis of neutral hydrogen gas is 
inclined by some 40 degrees to the axis of the inner disk, as 
revealed by the tilted ring model of Rogstad et al. (1976). 
More recently, Corbelli and Schneider (|1997p confirm that 
at a radius of 5 kpc (about 3 disk scalelengths) , the HI 
distribution in M33 shows a distinct change in inclination. 

Warps may be tidally generated when companions are 
obviously present; but in galaxies such as M33 which do 
not have any close companions, warping and gas infall may 



be inextricably linked (for a review, see section 7 in Binncy 
H992)) . 

In this paper, we focus our attention on the outer 
warped disk of M33, and in particular on its red stellar 
population. If the Hi-warp in M33 is induced by the in- 
fall of gas, then one interesting confirmation would be the 
presence of a red intermediate age (0.6 Gyr - 2 Gyr) popu- 
lation of carbon stars at the locale of the Hi-warp. Indeed, 
C-stars are expected to be associated with recent star for- 
mation. Tsalmantza et al. (2006) have found them towards 
the center in the LMC and associated with spiral arms in 
M31, the loci of star formation. In the SMC and in M33, 
their sample size was too small to trace their radial distri- 
bution, but Rowe et al. (2005) found them also relatively 
more abundant in the outer parts of M33. 

In an earlier study (Block et al. I2004[) . we presented 
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of red stars in the outer disk; the northern arc subtends 120° 
in azimuth angle and ~ 5' in width (see Figure [T]) . The 
northern arc is dominant although a very faint southern 
counterpart arc, forming a partial ring, can also be seen. 
The arcs lie at a radius of 2-3 disk scale lengths (in V, 
the disk scale length is 6 arcmin; Ferguson et al. [2006). 
Surprisingly, Fourier analysis of the deep 2MASS images 
showed that the dominant m — 2 peak did not correspond 
to any inner spiral arm morphology seen in optical pho- 
tographs, for which M33 is so famous. Rather, the m = 2 
peak corresponds to the giant outer red arms in the 2MASS 
images (see Figured]). 

Block et al. ; 2()() Q) noted that the very red colour of 
the arcs could not be due to extinction by dust grains. 
Freedman, Wilson and Madore (119911) derive a mean value 
for the total colour excess for the Cepheids in M33 of E(B- 
V) = 0.10 ± 0.09 mag, which includes both foreground 
(Milky Way) and internal M33 extinction. The K s -band 
extinction A(K S ) is approximately one-tenth that in the 
optical (Rieke & Lebofsky 1985]), so the dust extinction at 
K s is only ~ 0.03 mag. Extinction by dust grains cannot 
generate the very red colours found in the arcs in the outer 
disk of M33, where the extinction at K s would be even 
smaller. 

Our methodology differs from conventional ways of pho- 
tometrically identifying carbon stars using intermediate 
band filters (V, I, 77, 81) as presented by Rowe et al. 
(2005) for M33, using the Canada- France-Hawaii Telescope 
(CFHT). With such a method, one needs a large enough 
telescope to resolve individual stars. The (V, I, 77, 81) 
photometric method was pioneered by Cook, Aaronson and 
Norris (jl986|; another viable set of filters to identify car- 
bon stars are the Sloan filters. Demers & Battinelli (2005) 
used the three Sloan filters g, r, i, but again resolution of 
stars is crucial. Demers and Battinelli (2005) also used the 
CFHT with the Megacam camera to identify 361 new C- 
star candidates in M31. 

Some investigators (e.g. Tsalmantza et al. I2006P have 
focussed their attention on identifying carbon stars on the 
basis of the red J-K s tail (see Figures 1, 4, 6 and 9 in 
Tsalmantza et. al. 2006). The usefulness of such a method- 
ology is that near-infrared all sky surveys secured with 
moderate sized telescopes can potentially be exploited to 
identify carbon stars in spiral galaxies well outside our 
Local Group, where individual stars are not resolved. 

The arcs seen in Figure 1 are identified on the basis of 
their red J— K s colours. They lie around 15 arcmin in radius, 
between 12 and 18 arcmin, or 2-3 disk scale lengths, whereas 
the actual disk truncates at 5 scale lengths (Ferguson et al. 
2006). They are therefore like spiral arms in the outer disk, 
as seen also by Rowe et al. (2003|) . The latter authors have 
identified the C and M-type stars in the AGB population, 
and traced the Carbon star to M star ratio (C/M) as a 
function of radius. This ratio is not only an index of age 
but also of metallicity, increasing in metal-poor regions. 
Indeed, a C-star requires that the O-dominated surface be 
reversed to C-dominated, due to the dredged-up material 
from the star interior. This is easier to do when the surface 
is metal-poor. Rowe et al. (20Q5J) show that the C/M ratio 
increases with radius and then flattens beyond 20 arcmin 
(or 5 kpc), indicating a metallicity gradient covering most 
of the disk, including^ the outer snira.l arcs. 



of J— K s > 1.1. It was argued that while very old M gi- 
ants of solar abundance can indeed reach J— K s ~ 1 (see 
e.g. Figure 2 in Bessell and Brett |1998|) and even redder if 
they are super-metal-rich (see Frogel and Whitford [1987), 
stars with J— K s > 1 in low-metallicity regions cannot be 
M- giants but rather, very red carbon stars. There is a strong 
radial metallicity gradient in M33 (-0.16 dex/kpc in O/H, 
over 4-5kpc, see Beaulieu et al. 2006). Beaulieu et al. have 
confirmed this metallicity gradient, already found in HII 
regions, B-supergiants or planetary nebulae, by detecting 
beat Cepheids in M33. An excellent discussion of the metal- 
licity gradient in M33 may be found in Magrini, Corbelli 
and Galli (|2007j) . The outer regions of M33 are relatively 
metal-poor, and solar abundance is reached only in the very 
central domain of M33. 

In this paper, we examine the robustness of our photo- 
metric technique as a stepping stone to exploring the outer 
disks of more distant spiral galaxies, wherein individual car- 
bon stars may be present, but unresolved. To this end, we 
need follow-up spectroscopy. Here we discuss exploratory 
Keck I and Keck II spectroscopic observations of a few of 
the reddest and intrinsically brightest stars in the northern 
arc. All seven of our target candidates have very red J— K s 
colours and should be confirmed to be TP-AGB carbon 
stars if our methodology is correct. 

2. Observations 

As a prelude to conducting the Keck spectroscopy, we im- 
aged a section of the northern plume of M33 with the Hale 
5m reflector at Mount Palomar. We used the 2048x2048 
array near-infrared camera WIRC (Wilson et al. [2003), 
to resolve individual stars. The field of view is 8.5' x 8.5' 
with 0.25 arcsec pixels. The seeing FWHM is 0.8" in J, 
and 0.7 " in the K s band. The telescope was centered at 
01h34m28.1s, +30d54m00s (J2000). The total JHK S inte- 
gration time was ~ 9 minutes, reaching a limiting surface 
brightness at K s of 23.2 mag per pixel 2 , or 20.2 mag per 
arcsec 2 (4x4 pixels gives 1 square arsec). The point source 
photometry S/N — 10 limits are 19.0, 18.0, 16.9 mag in 
JHK S , respectively. 

Table [1] lists the magnitudes and colours of 7 candidate 
targets in the outer northern arc seen in Block et al. {2004) . 
Their positions are indicated by plus signs (+) in Figured 
In deriving the absolute magnitudes in Table [TJ we assume 
a distance modulus to M33 of 24.64'™, corresponding to a 
linear distance of 840 kpc (Freedman et al. I199ip . Noting 
that E(B-V) = 0.10 ± 0.09 mag (Freedman et al. ITWT1) . 
and that the K s -band extinction is approximately one-tenth 
that in the optical (Rieke and Lebofsky 11985ft . we use a 
dust extinction correction for the K s apparent magnitudes 
of 0.03 mag. 

Spectra of these carbon star candidates were obtained 
on 2004 Aug 17 UT with the Low Resolution Imaging 
Spectrograph (LRIS; Oke et al. 1995) attached to the Keck 
I telescope. The aperture mask slitlets had widths and 
lengths corresponding to 1'.'2 and 11'.'5, respectively. The 
600 groove mm -1 grating gave a resolution of ~ 8 A per 
projected slit width and a spectral range from 6825 to 9415 
A for slitlets near the centre of the mask. The stars were 
dithered along the slitlets between the 2 exposures, each 
6Dn seconds lone. The snectrnm from each slitlet, was ex- 
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Fig. 1. A partial ring of very red stars is seen in this J+H+K s image of M33, secured from a special set of 2MASS 
observations wherein the integration time was increased by a factor of six. The ring is seen to full advantage with a 
simple ellipsoid model subtracted. Black areas are those of excess emission while white areas correspond to areas of 
deficit emission. The northern plume or arc spans up to 5' in width and is located at a radius of 2 to 3 disk scale lengths 
(or 12 to 18 arcmin). An ellipse (of semi-major axis 16 arcmin) is drawn to pass through the arcs; North is up and and 
East to the left. 



calibration was produced by observing the spectrophoto- 
metric standard star G191B2B (Massey et al. I1998[) in one 
of the slitlets near the centre of the mask and transfer- 
ing the system response function found for this slitlet to 
the same wavelength intervals for the other slitlets. This 
approach neglects any differential vignetting or other field- 
dependent variations in response, but these are expected to 
be small for LRIS in any case. 

Also secured with the Keck II telescope were spectro- 
sconic observations of an area, fl.fi decree awav from the 



with very red J— K s colours in this outer low metallicity 
region. The data were obtained on 2004 July 18 UT using 
the Echellettc Spectrograph and Imager (ESI; Sheinis et al. 
2002) on the Keck II telescope. We used the echellete mode 
which allows a wavelength coverage from 3900 to 11060 A 
and yields a dispersion of ~ 56 km s _1 with the 1" slit. The 
spectra were obtained at parallactic angle near transit, and 
the total exposure was 1200 seconds for each object. 

The snectra. were reduced lisinc a. combination of .T. 
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Fig. 2. The positions of the seven stars listed in Table [T] are indicated by plus signs; the stars lie in the northern swath 
of carbon bearing stars identified photometrically in an earlier study (Block et al. I2004|) . Spectra of these stars (see 
Fig. [3|) were secured using the Keck I telescope in Hawaii. Also observed with its twin sister telescope, Keck II, are two 
stars which lie 0.5 degrees away from the centre of M33 in an area identified by a white circle. These outlying stars are 
representative of a very red family identified in near-infrared imaging with the Hale 5m reflector at Mount Palomar. In 
this DSS image, North is up and East is to the left. 



Table 1. Magnitudes and colors of our target stars in the 
northern red arc of M33 



Star 


J 


H 


K s 


J-K s 


J-H 


H-K s 


M Ks 


1 


17.96 


16.76 


15.91 


2.05 


1.20 


0.85 


-8.74 


2 


18.49 


17.09 


16.32 


2.16 


1.40 


0.77 


-8.33 


3 


19.61 


18.04 


16.66 


2.96 


1.57 


1.39 


-7.99 


4 


18.57 


17.26 


16.49 


2.08 


1.31 


0.77 


-8.16 


5 


18.82 


17.35 


16.54 


2.28 


1.47 


0.81 


-8.11 


6 


18.42 


17.24 


16.46 


1.97 


1.18 


0.78 


-8.19 


7 


18.37 


16.91 


15.82 


2.55 


1.46 


1.10 


-8.83 



ing spectrophotometric standards from Massey et al. (1998) 
also observed with the slit at the naralla.ctic anffle. 



Spitzer Space Telescope (Werner et al. 120041 Gehrz et al. 
2007) on 2005 January 21 as part of a Spitzer Guaranteed 
Time Observing Program (Program ID 5) conducted by 
Spitzer Science Working Group member R. D. Gehrz. 

The IRAC instrument (Fazio et al. 2004) is composed 
of four detectors that operate at 3.6/im (channel 1), 4.5/xm 
(channel 2), 5.8/Ltm (channel 3) and 8.0/im (channel 4). All 
four detector arrays are 256x256 pixels in size with mean 
pixel scales of 1.221, 1.213, 1.222, and 1.220"pixer 1 respec- 
tively. The IRAC filter band centers are at 3.548, 4.492, 
5.661 and 7.87 fim and the adopted zero magnitude fluxes 
are 280.9 Jy (channel 1), 179.7 Jy (channel 2), 115.0 Jy 
(channel 3) and 64.1 Jy (channel 4) as described in the 
IRAC Data Handbook V3.0 (see Gehrz et al. 2007). 

The M33 mannino' semience consisted of 438 frames ner 
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Fig. 3. Spectra of the seven stars seen in Figure [2j secured using the Keck I telescope. All C stars candidates are very 
red, with J— K s >2 mag; their individual colors are listed in Table [1] The coordinates of each star are indicated in each 
panel. Tick marks identify the prominent CN bands found in each star, shortward of 7000A and 8000A and longward of 
9000A. Each star is spectroscopically confirmed to belong to the thermally-pulsing asymptotic giant branch carbon star 
family. 



The raw Spitzer data were processed and flux calibrated software (Makovoz and Khan 2005]) • The artifact mitiga- 
with version 11.0.2 of the Spitzer Science Center (SSC) tion algorithm alleviates the effects of muxbleed, column 
pipeline. Post-BCD processing was carried out using an ar- pulldown/pullup, electronic banding, and bias variations 
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Detection, and Mosaicing. Background matching was per- 
formed by minimizing the pixel differences in overlapping 
areas with respect to a constant offset computed by the 
program. Cosmic rays and other outliers were detected and 
eliminated with the Outlier Detection module, fn the final 
step the images were reinterpolated to a pixel scale of ap- 
proximately 1.224"pixel _1 and mosaiced to create a final 
image spanning approximately 1.2° x 1.4°. 

3. Analysis of the Keck I and Keck II spectra 

Even for nearby M33, identification of carbon stars by 
spectroscopic means necessitates the use of the largest of 
groundbased telescopes. The first carbon star to be spec- 
troscopically confirmed in M33 dates back to the work of 
Mould and Aaronson (1986; Figure 7) who used the Hale 
5m reflector at Palomar. The individual V magnitudes of 
each target star in Table Q] are estimated, from their JHK S 
colours, to be fainter than 22 m , with one star (star 3 in 
Figure [3|) as faint as 25™. The spectra sample only 
the bright end of the K s -band luminosity function for car- 
bon stars, which ranges from Mk ~ -6 to -9. Much longer 
integration times would be needed for intrinsically fainter 
carbon stars. 

Carbon stars show a plethora of molecular spectral fea- 
tures, including the C2 Swan bands and the CN bands. 
The presence of the hugely dominant CN bands shortward 
of 7000A and 8000A and longward of 9000A are indicated 
by tick marks in Figure [3] and unambiguously reveal the 
C-star status for each of these seven stars. 
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4. Analysis of the ISM dust morphology in M33 

The carbon star arcs are clearly seen in JHK S images, but 
the global ISM dust morphology can also be effectively 
traced at 8/im. The reason is as follows: as carbon stars 
pulsate, their atmospheres become extended and matter 
may reach a temperature low enough for some elements to 
condense into very small particles ("star dust"). An out- 
flow develops as a result of radiation pressure and the star 
is progressively surrounded by an expanding circumstellar 
shell of dust and gas (Le Bertre et al. l2003|) . Emission from 
the shells of both carbon stars and also dust from O-rich 
shells is detected at 8[im (Buchanan et al. l2006|l . Buchanan 
and co-workers note that "most of the C-rich stars have 
spectra that are dominated by warm dust" and these au- 
thors show that stars such as MSXLMC587 (O-rich) and 
MSXLMC1488 (C-rich) both show emission in the IRAC 8 
fim window (see their Figure 2, for example). The trend is 
for the flux to climb from longer (~ 38^m) to the shorter 
wavelength of 8/zm. As an interesting aside, one does not 
need high ultraviolet excitation to produce infrared emis- 
sion features: stars as cool as 3000 K can generate PAH 
features, for example, as discussed by Li and Draine (2002). 

Carbon stars are one of the most important contribu- 
tors to the replenishment of dust in the interstellar medium. 
Our IRAC image of M33 at wavelength 8/im offers un- 
precedented insight into the global ISM dust morphology. 
Deprojection of our IRAC images of M33 were facilitated 
by adopting log R25 = 0.23(RC3) where R is the ratio of 
the ma.ior tn minor axis, and a. nosition ancle nf 93 decrees 



Fig. 4. Spectra with the Keck II telescope, of two stars 
0.5 degree (~5 scalelengths) away from the centre of 
M33 (the field circled in Figure \2§ , centered at approxi- 
mately RA(2000)=lh34m21s and DEC(2000)=31°08'45"7. 
The dominant CN absorption features are clearly evident 
and betray the signature of TP-AGB carbon stars. 



the Rayleigh- Jeans tail) from the 8 /j,m image, we subtract 
a scaled version of the IRAC 3.6 fj,m map. As in Block et 
al. (2006), we generate a 'non-stellar' pure dust map from 
8/xto - k 3.6/iTO, where n (a constant) assumes the value of 
0.25. This subtraction removes the contribution from stel- 
lar photospheres and leaves only the emission from dust 
grains, which trace the ISM morphology. The major uncer- 
tainties entering into the precise value of k are, firstly, the 
extended aperture correction, and secondly, the exact mid- 
IR colour of the stellar populations which one is modeling 
by scaling the 3.6 /im emission. The second correction in- 
troduces uncertainties of order 0.05 mag or perhaps larger 
(Ashby, private communication) . A slightly smaller value of 
k (0.232 instead of 0.25) has been suggested by Helou and 
co-workers p004[) . but our use of a slightly higher photo- 
spheric contribution (at a level of 2 percent) is well within 
the correction uncertainties; the choice of k is only accurate 
at the 10-15 percent (and not one or two percent) level. 

A complimentary method could, of course, be to focus 
our attention on photosphcric emission, and not ISM dust 
emission. At the shorter IRAC wavelength of 3.6/im, car- 
hon star nhotosnheres would he detected, hut sn would the 
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Radius (arcmin) 



/ i{u,e)dude 

-7T J — OO 

In practice, the integrals in u = In r are calculated from 
a minimum radius (selected to exclude the bulge where 
there is no information of the arms) to a maximum ra- 
dius (which extends to the outer limits of the arms in our 
images). 

The inverse Fourier transform can be written as 



M= 2 




10 20 
Radius (arcmin) 

M= 2 




Radius (arcmin) 

Fig. 5. The topmost panel shows a dust emissivity 8/ito 
radial profile. In the middle panel is shown the amplitude 
of the Fourier m=2 mode as a function of radius. The lo- 
cal peak at r — 14' coincides with the inner radial limit 
of the arcs. The outer spiral arms or "arcs" correspond 
to the higher m—2 amplitude between 14 and 18 arcmin. 
The phase of the m=2 component is presented in the lower 
panel. We see a change of phase in the dust emissivity distri- 
bution at 18 arcmin, at the outer boundary of the outer spi- 
rals, precisely the same radius where Corbelli and Schneider 
ftW7\i find a warp in HI. 



(Gehrz et al., in preparation). In this paper, we confine our 
attention to the dust grain morphology, as we did for the 
Andromeda Spiral in Block et al. ((2006). 

The Fourier method has been extensively discussed in a 
number of papers (e.g., Considere and Athanassoula:1982; 
Puerari and Dottori 119921 among others) . In the Fourier 
method, an image is decomposed into a basis of logarithmic 
spirals of the form r=r Q exp(— ^8). The Fourier coefficients 
A{p 1 to) can be written as 



1 l"^ 1T f^ 00 
MP: m ) = 77 / / I ( u > °) ex P 

U J — 7T J — OO 



-i(m9 + pu)]dud0 



Here, u = In r, r and 6 are the polar coordinates, m 
represents the number of the arms, p is related to the pitch 

an (rip P of tho snira.l hv P=atan( — m. Ini\ and T(n f)\ is 



where 



S m {u) 



S(u,6) 



D 

!«47r 2 



S m (u)e 



iinO 



G m (p)A(p,m)e* pu dp 



and G m (p) is a high frequency filter used to smooth the 
A{p, to) spectra at the interval ends (see Puerari and 
Dottori 1992), and it has the form 



G m (p) = exp 



1 fp-p 
2 



25 



where p™ ax is the value of p for which the amplitude of the 
Fourier coefficients for a given to is maximum. The chosen 
interval ends (p+=+50 and p-=— 50), as well as the step- 
size dp=0.2b, are suitable for the analysis of galactic spiral 
arms. 

Firstly, the Fourier m = 2 mode in the 8/ito image also 
does not correspond to the inner spiral arm morphology, 
but rather shows a local peak at a radius of 14', just long- 
ward of 2 scalelengths. This is precisely where the outer red 
arcs were detected by Block et al. p004p . 

Another intriguing fact is that there is a distinct change 
of phase in the ISM dust morphology at 3 disk scale-lengths 
(18'), which is precisely the region of the outer bound- 
aries of the "arcs" or outer spiral arms, where Corbelli and 
Schneider (1997) find the prominent Hi-warp. 

We believe it fully consistent to propose that the very 
red, and relatively metal-poor stars have recently been 
formed by gas infall which is inextricably tied to its strong 
Hi-warp. It seems highly plausible that fresh low-metallicity 
gas is being fed to the host galaxy M33. The outer disk from 
which mass is accreted is inclined to the inner disk of M33, 
and has a different angular momentum, as revealed by the 
tilted disk model. 

Our conclusions are compatible with the work of Rowe 
et al. (|2005p who published a spatial map of the C /M-star 
ratio, that reveals conspicuous maxima at exactly the posi- 
tions of what we call "arcs" , or outer spiral arms. They in- 
terpret the radial distribution of the C /M ratio uniquely as 
a metallicity indicator, and conclude that there is a strong 
metallicity gradient up to 5kpc (20 arcmin) and then the 
gradient flattens out. Their interpretation is in terms of the 
viscous disk model, where the outer shear of the rotation 
curve produces mixing in the outer parts, able to wash out 
abundance gradients. 

The viscous disk model certainly applies to M33 to some 
extent, and is certa.inlv a. needed incredipnt of tine external 
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The viscous accretion disk model has been frequently in- 
voked to build exponential galaxy disks, either with viscous 
torques to exchange angular momentum, or with gravita- 
tional torques due to spiral arms or bars. In the latter case 
the equivalent viscosity is called "gravitational viscosity" 
(e.g. Lin & Pringle fl987|) . If the viscous and star forma- 
tion time scales are of the same order, then the resulting 
stellar disk has an exponential distribution independent of 
the disk rotation law and of the assumed viscosity prescrip- 
tion. The expected metallicity gradient is also exponential, 
as computed by Tsujimoto et al. (flD95) . As for M33, it is 
highly possible that the "viscous disk" model applies until 
a radius of 5kpc (20 arcmin), i.e. just outside the arcs, and 
is the mechanism with which the accreted gas coming from 
the warped component is consumed into stars and disk for- 
mation. The optical disk is exponential up to this radius. 



5. Conclusion 

We have confirmed, through Keck spectroscopy of individ- 
ual stars, the presence of an enhanced abundance of carbon 
stars in the outer spiral arms in M33, in particular in the 
two "arcs" identified previously from their red colour (Block 
et al. 2004). Our conclusions are in very good agreement 
with the recent work of Rowe et al. (2005), who have pro- 
duced a spatial map of the C /M-star ratio in M33 showing 
the conspicuous outer arms. We also see the outer spiral 
arms in the Fourier m=2 component of the Spitzer 8fim 
dust emission image. These arms correspond to the radius 
where the HI-21cm warp starts, and we propose an inter- 
pretation in terms of recent star formation, fueled from the 
gas accreted from this outer reservoir. A very similar con- 
clusion, in terms of M33 accreting gas, has recently been 
reached by Magrini, Corbelli and Galli (200T|), on the ba- 
sis of O/H, S/H and [Fe/H] abundances in the Triangulum 
Spiral Galaxy M33. 

The HI gas reservoirs in the outer parts of galaxies are 
observed outside nearly all spiral disks (Sancisi ll983p . Their 
conspicuous warped morphology, even in the absence of any 
perturbation or companion, implies that this gas is being 
almost continuously accreted, with a different angular mo- 
mentum than that of the inner disk (e.g. Binney 1992). 
Keres et al. (|2005|) discuss cold gas accretion along fila- 
ments in the cosmic web. Accreting systems in gas need 
not show any signs of accretion in stars, such as the pres- 
ence of tidal tails, stellar loops or close companions. 

The implications of the presence of carbon stars in the 
outer disk of M33 immediately beckons the question of a 
possible ubiquity of such an intermediate age population 
in the outer domains of other spiral disks out to the Virgo 
cluster and beyond. 

At hand is the potential to apply our photometric 
technique of identifying thermally-pulsing asymptotic gi- 
ant branch carbon stars to more distant spirals on the ba- 
sis of their red J-K s colours in low metallicity domains, an 
approach adopted by other teams such as Tsalmantza et 
al. ([2006). It is a sobering thought that it takes a ground- 
based class 8— 10m telescope to secure individual carbon 
star spectra in the outer disk of M33, our second closest 
spiral, in realistic amounts of observing time. Hence the 

nrpent need to identify unresolved sets of carbon stars in 



class telescopes, and then imaging them through Sloan or 
(V, I, 77, 81) filters. 

Carbon stars may make an important contribution to 
the IR luminosity of high-redshift galaxies. Using the Large 
Magellanic Cloud as a guide, carbon stars are produced 
in large numbers between ages of about 0.6 and 2 Gyr. 
Therefore galaxies that undergo a burst of star formation 
will have their IR light boosted 0.6 Gyr later, and this ex- 
tra light will die away after about 2 Gyr. We know that 
star formation in the early universe (Bouwens et al. [2003 ) 
was already proceeding at redshift of at least z ~ 6. There 
may be an epoch starting about 0.6 Gyr after the onset of 
star formation in the universe (i.e. at redshifts z ~ 4) char- 
acterized by large numbers of carbon stars. The dominant 
output from the carbon stars will be redshifted into the 
mid-infrared where these stars could double the observed 
extragalactic flux. Maraston (|2005|) indeed proposes the use 
of carbon stars as an age indicator for high-redshift stellar 
populations. Instruments like the Mid-Infrared Instrument 
on board the JWST are needed to image such galaxies in 
their (rest-frame) 2.2 micron band. 
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